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1. Introduction 

Heavy quark physics presents the lattice with an invaluable opportunity to provide both theo- 
retical insight and the key to some precise results from experiment. The use of effective field theory 
techniques such as NRQCD [[Tj] allow heavy quarks to be described on the lattice avoiding the large 
discretisation errors of relativistic formulations. 

The basic idea of these techniques is that simulations should be performed on lattices coarse 
enough that the cutoff n/a separates the heavy quark mass 1 scale M from the dynamical scales 
Mv, Mv 2 , where v is the heavy quark velocity inside the bound state. This allows a non-relativistic 
formulation of QCD to simulate only the dynamical scales important to heavy mesons, with the 
physics at the heavy quark mass scale (and above) being encoded into matching coefficients. 

The other key component of the calculations presented here is the use of 2 + 1 flavour dynami- 
cal configurations provided by the MILC collaboration [Q]. These configurations employ the a 2 im- 
proved "AsqTad" improved staggered quark action and a highly improved gluon action to describe 
the vacuum. Our calculations look at the coarsest ensemble generated which is comprised of around 
600 16 3 x 48 lattices with spacing a ~ 0. 17 fm, intended to complement the coarse (a ~ 0. 12 fm) 
and fine (a ~ 0.09 fm) ensembles on which the majority of results have been extracted. Part of the 
motivation for our calculations was to determine the feasibility of the extra-coarse ensemble for 
extraction of physical results in the presence of possible large discretisation errors. 

In section ||| we outline the formulation used for our calculations, then in section || we look at 
the T and y spectra. In section |] we calculate the mass of the B c meson using NRQCD for both the 
b and c quarks. Finally, in section |5] we look at the B s and Bj mesons by combining an NRQCD b 
quark with an improved staggered light quark. 

2. Formulation and Technical Details 

A potential model calculation using a logarithmic potential |^] can be used to extract the 
typical velocity of the heavy quarks in their bound states. Doing so gives v\ ~ 0.1 for bb, v 2 = 0.25 
for cc and v\ ~ 0.04, v 2 ~ 0.4 for the B c meson. These values justify a non-relativistic description, 
with the possible exception of the c quark of the B c system. This work makes use of the non 
relativistic formulation of QCD (NRQCD) developed in [j|]. In this formulation, the operators of 
the heavy quark hamiltonion are classified according to a power counting in the velocity of the 
heavy quark in the bound states; v. We use the action given in[JT|] which is correct through €?(v 4 ). 

We fit the mesonic corellators using the constrained curve fitting method of [Bh. This method 
employs Bayesian techniques which utilise prior information to guide the fit, allowing all of the 
excited states present in the initial few time slices of correlator to be fitted. This improves the 
statistics of the fits while reducing the sensitivity to the number of exponentials used. 

3. The T and y/ heavy quarkonium systems 

The majority of calculations using NRQCD describe b quark systems, see for example [|0. In 
the quantity ciMq, Mq for the bb is a relatively large number {Mq » 1), so that even with present 

'Throughout this work, the term heavy quark refers to either the b or c quark. 



217/2 



NRQCD results on the MILC extra coarse ensemble 



Ian Allison 



day lattices aM q is quite large. On the same lattices, the cc system requires an Mq of less that 
1.0 and inverse powers of Mq in the NRQCD evolution equation then cause numerical instabil- 
ity, preventing results from being extracted. The advantage of the extra coarse configurations is 
that Mq = 1.0 produces mesons whose kinetic mass is close to that of the particle while the 
evolution equation remains stable. 

In NRQCD, the zero of energy is shifted, so that the energy of a meson with momentum p is 



E(p) =E(0) + Vp 2 + M 2 - M, 



(3.1) 



which allows us to extract a physical mass using states at non-zero spatial momentum. This kinetic 
mass is commonly used to tune the input quark mass. From figures [| and |2| it is clear that the b 
quark mass (Mb = 4.0) is accurately tuned, while the c quark mass (M c = 1.0) is about 12% too 
low. 
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Figure 1: Kinetic Mass of the T 



Figure 2: Kinetic Mass of the J/\jf 



Rearranging expression (|3jJ) it is possible to extract the effective speed of light; c = ((AE) + 



2MAE) I p 2 , where M is the kinetic mass determined at some fixed p (taken to be p = 1 in this work). 
Of course, this quantity should be equal to 1, but on the lattice discretisation effects can ruin this. 
We found c 2 = 1 to within statistical errors of 1 % and 4% for the bb and cc systems respectively. 
This compares very favourably with relativistic formulations without high levels of improvement, 
where deviations of up to 10% are not uncommon 

The systematic way in which NRQCD was developed makes it possible to estimate the effect 
of higher order terms left out of the simulation hamiltonian. To do this the expectation values of 
these operators are calculated using potential model estimates of (p 2 ) and (p 4 ). These corrections 
fall naturally into three categories. 

Relativistic corrections: The main relativistic correction comes from the first term left out of 
the relativistic expansion of the energy momentum relation 2 , (p) 6 /8Mg. 

Radiative corrections: The effective field theory coefficients c,- of the terms appearing in the 

NRQCD hamiltonian take the form c,(Mg,a s ) = l + c ( . (Mo)a s -\ , where the c ; - are calculable 

in perturbation theory but are expected to be £?(1). We estimate the effect of radiative corrections 
for the kinetic term; a v (p 4 ) /4(Mg) 3 , and for the Darwin term 3 ; 47ra 2 i// 2 (0) 2 /(3Mg) 2 . 

2 The correction terms all have an extra factor of 2 included because of the two quarks in the meson. 
3 This expression comes from estimating the A.E term in a potential model. 
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Discretisation corrections: The terms in the heavy quark hamiltonian which correct for dis- 
cretisation errors in the lattice spatial derivative and temporal derivative have corrections of the 
form a s a(p 4 )/Sn(M Q ) 2 and a s a 2 (pf)/ 12 (M Q ) respectively. 

In this work we are interested only in energy splittings, so that adding errors from these sources 
for each state in quadrature is an over estimate. Instead, we calculate the above corrections using 
potential model estimates of the difference in the expectation values of the various splittings. The 
values calculated are shown in table [I]. It can be seen from this table that the errors are substantially 
less for the radial splittings than for the orbital splittings. This is just a statement of the similarity 
of the dynamics of the radial states leads to a more efficient cancellation of higher order effects. 



Correction 


T(25- IS) 


y(2S- IS) 


t(ip- is) 


\I/{IP-IS) 


a,V/4(M e ) 3 


0.55 


1.07 


2.19 


1.73 


Sp 6 /(M Q ) 5 


0.30 


1.64 


0.92 


1.47 


47ra s 2 y/- 2 (0)/3M 2 


0.62 


0.90 


1.79 


2.82 


Total ReL/Rad. 


0.88 


2.16 


2.97 


3.62 


a s a 2 8pj/l2M Q 


2.2 


0.61 


8.76 


0.99 


a s a8p A j%M 2 n 


0.61 


0.23 


2.45 


0.37 


4na s a\i/ 2 (0)/l5 


2.47 


0.51 


7.17 


1.60 


Total Disc. 


3.36 


0.83 


11.58 


1.92 


Total 


3.47 


2.31 


11.96 


4.09 



Table 1: Systematic corrections to heavy quarkonium splittings as a percentage of the experimental splitting 

From these determinations it is clear that the 25 — IS splitting of the T system presents the best 
opportunity for extracting the lattice spacing for bb states, it gives a~ l = 1.12 1(62) (31) where the 
first error is statistical and the second comes from the above systematic considerations. Increased 
statistical precision means that IP — IS) splitting the best choice for setting a from the cc system, 
giving aT x = 1.111 (29) (23). The agreement of these two determinations is unsurprising when 
using dynamical gauge configurations. 

In figures |3| and || we plot the spectra of bb and cc systems respectively. The agreement of 
the lattice points with the experimental values (shown as lines) is good in the T system. In the \ff 
system, there is a noticeable discrepancy in the hyperfine splitting of 44 MeV. This is too large to 
be fully explained by the systematic corrections of table |l| but is at least in part due to the fact that 
the c quark mass was tuned too low. 

4. The B c Mass 

Last year we presented a calculation of the mass of the B c meson using NRQCD b quarks and 
Fermilab method c quarks [|7|]. The use of dynamical configurations in this calculation allowed a 
level of precision which, when combined with the experimental status of the B c mass at that point, 
allowed us to call our calculation a prediction. We found Mg c = 6.304(4)(11)(+q 8 ) [B], where, in 
order, these errors were due to statistics and chiral extrapolation, uncertainty in setting the input 
quark masses, and an estimate of higher order corrections to Fermilab quark action. 
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'S (0'*) 3 S, (1") 1 P, (1*5 3 PSpinAvg. 3 P(0,1,2") 




S (0'*) 3 S,(r) 'P, (1*") 3 PSpinAvg. 3 P (0,1,2" 



Figure 3: The T spectrum. X P( IP - 15) has been Figure 4: The y spectrum. The 15) has been 

used to set the lattice spacing and the absolute scale used to set the lattice spacing and the absolute scale 
has been set by the 3 5i has been set at the spin average of the 15 states 



Soon after our result was published, new determinations of the B c mass were made at the 
Tevatron [Q]. They quoted a preliminary result of Mb c = 6.2879(48)(11) where the errors were 
statistical and systematic respectively. 

Our method relies on forming the difference, 

m Bc - 1/2 (m T + m J/¥ric ) = M Bc - 1/2 (M T + M J/¥Vc ) (4.1) 

where the quantities on the left are lattice energies, and those on the right are physical masses. The 
states labelled yr\ c are the spin average of the \ff and r\ c . Repeating our earlier calculation on the 
extra coarse ensemble but using NRQCD to describe both valence quarks, we find 6.274(1)(25), 
where the first error is statistical, and the second is calculated in the same way as those in table |]. 
Encouragingly, this value is in agreement with both the experimental value and our previous result. 



5. Heavy light results 

Finally, we also looked at the B s , and Bj systems using and NRQCD b quark and an improved 
staggered light quark. This calculation is complicated by the fact that the fitted correlators include 
an oscillatory contribution from the parity partner of the desired meson. This is a consequence of 
using staggered light quarks. 

The results we present here are preliminary, and the reader is referred to [ JT0[ ] for a more com- 
plete treatment. In figure [5] we plot the quantity 2mg s — «ir ■ Terms explicit in the quark mass cancel 
in this difference, exposing the differences in the binding energies of the systems. Our determi- 
nation of 2mg s — my seems to be systematically too large, but we present here only a statistical 



error and delay a systematic determination of the errors until [JlOp , where further results will also 
be presented. It should however be noted that the sea s quark mass included was estimated to be to 
high by approximately 30%. 

As yet our results for the B s and Bj do not support any conclusions, but calculations increasing 
the number of light quark masses and simulating the D s and D with NRQCD c quarks are currently 
being analysed. Given the success of our heavy quarkonium calculations, there is a good prospect 
that these calculations will yield useful results. 
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Figure 5: 2ms s — m v . The boxes are our results 
and the star is experimental value placed at the 
physical s quark mass (am, = 0.062) for this en- 
semble. 
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Figure 6: The Kinetic mass of B s in physical 
units vs. spatial momenta. The errors shown 
are purely statistical. Experimentally Mb s = 
5.3696(24) GeV 
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